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ABSTRACT 

We report on two optical candidates for the counterpart to an X-ray source in the SMC, 
IWGA J0053.8-7226, identified as a serendipitous X-ray source from the ROSAT 
PSPC archive, and also observed by the Einstein IPC. Its X-ray properties, namely 
the hard X-ray spectrum, flux variability and column density indicate a hard, transient 
source, with a luminosity of ^3.8 xlO"^^ergs sec~^. XTE and ASCA observations have 
confirmed the source to be an X-ray pulsar, with a 46 s spin period. Our optical 
observations reveal two possible candidates within the error circle. Both exhibit strong 
Ha and weaker Hf] emission. The optical colours indicate that both objects are Be-type 
stars. The Be nature of the stars implies that the counterpart is most likely a Be/X- 
ray binary system. Subsequent IR photometry (JHK) of one of the objects shows the 
source varies by at least 0.5 mag, while the (J-K) measured nearly simultaneously with 
the UBVRI aivd spectroscopic observations indicate an IR-excess of ^0.3 magnitudes. 

Key words: stars: emission-line. Be - star: binaries - infrared: stars - X-rays: stars - 
stars: pulsars 



1 INTRODUCTION 

Accretion onto a compact object invariably results in X-ray 
emission, with the luminosity determined by such important 
parameters as the size and mass of the accreting object, 
and the mass accretion rate. Binary stars represent the best 
laboratories for the study of such accretion processes, where 
the accreting star can be a white dwarf, neutron star or 
a black hole. Whereas known white dwarfs in cataclysmic 
variables are generally low luminosity X-ray sources, and 
relatively nearby (most <1 kpc), the neutron star sources are 
often extremely luminous, and detectable at extragalactic 
distances. 

Most (~80%) of the higher mass systems (Ifigh Mass 
X-ray Binaries, or HMXBs), which consist of early spectral 
type (O or B) stars losing mass to their neutron star com- 
panions, are Be/X-ray or Supergiant binaries. The orbits of 
these binary systems are generally wide and eccentric, with 
periods in the range 0.7/1.4d for RX J0050.7-7316 (Coe & 
Orosz, 2000) to 262d for SAX J2239.3-h6116 (in't Zand et 
al, 2000). Usually the only optical signature of these stars 
is Balmer emission, and in that respect alone it is often dif- 



ficult to ascertain the X-ray nature of these systems. With 
Be stars being relatively common in the Galactic plane, it 
has only been possible to confirm Be counterparts to X-ray 
sources using imaging X-ray instruments, as provided by 
such sateUites as Emstem (HEAO-2) and ROSAT. 

Much progress in understanding the physics of these 
systems has resulted from multiwaveband campaigns, span- 
ning the IR-optical-UV-X-ray-7-ray domains. Observa- 
tions of the Be star in the optical and IR sheds light on the 
physical conditions under which the neutron star accretes. In 
addition, with contemporaneous X-ray observations, which 
provide direct information on the accretion process, it is 
possible to investigate the correlation between the mass loss 
rate from the Be star and the accretion rate. Thus it is possi- 
ble, through long-term multi-wavelength programs, to build 
a more complete picture of the whole accretion process, and 
to study it as a function of time. 

In this paper we report on two possible counterparts to 
a hard X-ray transient source in the SMC, IWGA J0053.8- 
7226, which is most likely a new member of the Be/X-ray 
binary subclass of HMXBs. The .ROSAT' position coincides 
with a previously know hard spectrum source from HEAO-2 
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[Einstein Observatory) IPC pointings. Our optical/IR obser- 
vations of one of the candidates reveal IR variability of ^^0.5 
mag and persistent Balmer emission, with a redshift consis- 
tent with SMC membership. The UBVRI colours and optical 
spectra of both candidates indicate an early B-spectral type. 

In 1997 November, following our observations, an X- 
ray source identified as XTE J0053-724 went into outburst 
and revealed a 46 s X-ray pulsation (Corbet et al. 1998). 
Though originally, and incorrectly, identified as SMC X-3, 
subsequent ASCA observations in 1997 December revealed 
that the pulsating source was actually 8' away from SMC 
X-3, establishing it to be a new X-ray pulsar. Its new ASCA- 
determined position is coincident with the previously cata- 
logued X-ray source IWGA J0053.8-7226 and the data we 
present here confirms this association and provides convinc- 
ing evidence that it is a Be/neutron star binary. 
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2 ROSAT OBSERVATIONS 

The WGA catalogue (White, Giommi, Angelini 1994) was 
generated using all pointed ROSAT PSPC data available 
in the archive up to 1994 November and its first revision 
includes data up to 1995 March. The catalogue was cre- 
ated using the revl ROSAT data files, where data taken 
several months apart belonging to the same sequence were 
combined in one file. A variability test for each detected 
source was included as part of the processing. The method 
used compares the time arrival distribution of the photons 
collected in each pixel with the corresponding distribution 
of the entire image using a Kolmogorov-Smirnov (KS) test 
(Giommi, Angelini, White 1995). The transient nature of 
IWGA J0053. 8-7226 was discovered using this method, in 
the sequence 600195, which included data from 1991 Octo- 
ber and 1992 April. The source was detected in 1992 April 
at coordinates R.A.=00'' 53™ 53.8' and Dec=-72° 26'35" 
(J2000). Positional errors from the WGA positions vary from 
13 — 50", depending how far the source is off-axis. 

After discovering the transient behaviour, the ROSAT 
public archive was searched for additional PSPC and HRI 
public data containing the source position. Table 1 gives the 
observation log Q together with the nominal exposure, the 
off-axis position of the source in the HRI or PSPC Field of 
View measured from the nominal pointing, and the count 
rate. 

A count rate (or upper limit) is given in Table 1 for 
most of the observations with the following caveats. For a 
large ofi^-axis angle, the point spread function (psf) is very 
large: at 50 arcmins the 50% power radius is ~ 3 arcmins 
(see ROSAT Handbook). That far off-axis, the source posi- 
tion is typically at the detector boundary and the optimized 
radius falls outside the detector image, resulting in an un- 
derestimate of the countrate. IWGA J0053. 8-7226 is at ~ 5 
arcmin from another transient RX J0054.9-7210 (Kahabka 
and Pietsch 1996), which is also an Emstem IPC source: No. 
9 in Bruhweiler et al (1987) and No. 35 in Wang and Wu 
(1992). At such large off-axis angles, if this second source 

* We used the rev 2 data processed with SASS7, which included 
among other changes, a better attitude reconstruction, and the 
spUtting of the sequence into separate files when the guide stars 
change (see also: Pisarski 1995). 



Figure 1. The background subtracted lightcurve is obtained for 
the PSPC observations 600195a01 and 600196a01 for the 1992 
April observation. The data are corrected for PSF and vignetting 
effects. Each bin corresponds to 1000 seconds. 



was active, the optimized radius to evaluate the count rate 
for IWGA J0053. 8-7226 is contaminated by emission from 
RXJ0054.9-7210. For these observations no count rate or up- 
per limits are given in Table 1, although a column indicates 
whether IWGA J0053.8-7226 was detected ('on' or 'off'). 

The source is highly variable and it was detected in 
1992 April (PSPC), 1993 October (PSPC), and 1996 June 
(HRI). Although the coverage is not complete, the outbursts 
for the 1992 April and 1993 October observations last more 
than a month. Figure 1 shows the background subtracted 
lightcurve for the 1992 April outburst combining data from 
the sequence 600195 and 600196, where only the declining 
phase is seen. The PSPC spectrum obtained for the 1992 
April observation (only when the source is at ~ 20 arcmin 
off-axis, because of the response uncertainties at large off- 
axis angles) is shown in Figure 2. We fit the spectrum using 
a power law and absorption. The best fit parameter values 
are a spectral index of 1.6 with an equivalent hydrogen ab- 
sorption of 3 X 10^^ cm~^ (shown in Figure 2), but these two 
parameters can not be simultaneously constrained. 

A model with a fixed absorption at the galactic value 
in the SMC direction ( 6 x 10^" cm~^) and a power 
law, overestimates the low energy part of the spectrum. If 
IWGA J0053.8-7226 is a Be star, as argued later, it is rea- 
sonable to assume that the outburst spectrum should be 
similar to that of IWGA J0051.8-7231, another Be tran- 
sient in the SMC (Israel et al. 1997). Fixing the spectral 
index to 1.1, as for IWGA J0051.8-7231, we obtained an 
N_r/=(2.3 ± 1.1) X 10^^ cm~'^. This value seems consistent 
with the contour map in Kahabka & Pietsch (1996) where 
IWGA J0053.8-7226 is close to source number 103 in their 
paper. This gives a 0.1-2.0 keV unabsorbed X-ray flux of 
7.5 X 10"^^ erg cm^^ s~^ which corresponds to a luminosity 
of ^ 3.8 X 10''^ erg s~^ for a distance of 65 kpc. 

The 1996 June HRI observation allowed us to derive a 
better position because of the higher HRI angular resolution 
compared to the PSPC, and because the source location is at 
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data and folded model 



IWGA J0053, 8-7226 



N 



Figure 2. X-ray spectrum, power law fit and residuals for the 
600195a01 observation oflWGA J0053.8-7226. 



a smaller ofT-axis angle (~ 3 arcmins) compare to the PSPC 
observations. At equinox 2000 this is R.A.=00'' 53"" 55.6^ 
Dec=— 72° 26'44" with a conservative error of 15" (see also 
HRI report, Harris (1996)). 

The position obtained from the PSPC No. 3 ob- 
servation using the rev2 data is R.A.= OO*" 53™ 56.0° 
Dec=-72° 26' 34.8" (±25"). The IWGA J0053.9-7226 po- 
sition is consistent with the Emstevn IPC source No. 8 
in Table 2 of Bruhweiler et al (1987) and source No. 34 
in Wang and Wu (1992). The coordinates of the Einstein 
source from Wang and Wu (1992) are RA=00'' 53"" 54.2= 
Dec=-72° 26' 27" in equinox 2000, with an error ±40". 

The source was noted by Bruhweiler et al. to have the 
hardest spectrum amongst the 35 listed, and assumed to be 
a stellar source. From an IPC count rate of 0.0088 cts s~^, 
Bruhweiler et al. determined an X-ray luminosity of 1.8 
xlO^'^ergs s-i and hardness ratio HR ={H- S)/{H ± S) 
= (0.8 to 3.5 keV - 0.2 to 0.8 keV) / (0.8 to 3.5 keV ± 0.2 
to 0.8 keV) = 0.76, consistent with the hard spectrum we 
report here. 



3 OPTICAL IDENTIFICATION 

The field of IWGA J0053.8-7226 is shown in Figure 3. There 
are 4 stars lying within the X-ray error circle labelled as 
Stars A, B, C and F. Star A is, in fact, a double and both 
components were subsequently observed on 1995 September 
22/23, using the SAAO 1.9-m telescope with the Cassegrain 
spectrograph. Spectra in the interval 3300-7500A, at a res- 
olution of ~7A, were obtained using the Reticon Photon 
Counting System (RPCS). Because the two components of 
Star A were separated by only a few arcsec in the E-W di- 
rection, it was impossible to get separate spectra of the two. 
The combined spectrum was clearly very blue, and with a 
strong Ha emission line (E.W. = -13A). Although the RPCS 
is only a ID detector, it was possible to determine which 
of the two stars was responsible for the Balmer emission. 
This was done by driving the stellar image off to one ex- 
treme of the slit decker, thereby excluding most of the light 
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Figure 

IWGA 



3. Finding 
J0053.8-7226 



chart 
field 



(~2 X 2 arcmin) for the 
from the Digital Sky Survey 
(http://archive.stsci.edu/dsf) showing all the possible counter- 

The HRI 90% confidence error cir- 
the centre of which is R.A. = 
(Equinox 2000). 



parts referred to in the text 
cle (radius = 15") is marked 
OOHsa^SS.e^Dec = -72°26'44' 



of the companion. The exercise was then repeated for the 
other star. It was therefore possible to obtain two spectra, 
whose major light contribution was one or other of the two 
stars. From these observations we were able to show that 
the fainter, easternmost, star was the Ho emitter. This has 
been subsequently confirmed from both optical photometry 
and CCD spectroscopy, reported later in this paper. 

CCD photometry of the field of IWGA J0053.8-7226 
was made using the SAAO 1.0-m telescope on 1995 Novem- 
ber 26/27. The observations were made using the Tek8 CCD 
camera, employing UBVRcIc and Ha filters, at a plate 
scale of 0.35 arcsec pixel"^. All observations were obtained 
in photometric conditions, and in good seeing (FWHM 
-1"), with exposure times of 600 s (17), 300 s (5), 200 s 
{V,Rc & Ic) and 1000 s (Ha). 

All of the CCD frames were reduced in a similar man- 
ner, with bias removal, flat-fielding, and subsequent reduc- 
tions using DoPHOT (Mateo and Schechter 1989). Obser- 
vations of E-region standard stars (Menzies et al. 1989) 
were used to determine aperture magnitude zero points. 
Both the PSF fit and aperture magnitudes were determined 
for all stars in the frames, and transformed to the natural 
UBVRcIc system. Colour-magnitude and colour-colour di- 
agrams were made in order to search for stars of unusual 
colour. In the past the Ha filter has been used successfully 
in conjunction with the 7?-band filter, to search for Ha emis- 
sion line objects. No stars on the frame appeared to show 
significant deviations suggestive of very strong Ha emission. 

Our CCD observations resolve the two components of 
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Table 1. Table 1: Observation log of IWGA J0053.8-7226 



No. 


Instrument 


Year 


Day /Month 


Exp.i 


Seq. 


OfT-axis 
(arcmin) 


On/Off 


Count rate ^ 
(xlO-2) 


1 


PSPC 


1991 


8/Oct 3/Nov 


16644 


600195AOO 


21.9 


off 


< 0.536 


2 


PSPC 




9/Oct 2/Nov 


1303 


600196AOO 


49.3 


off^ 


- 


3 


PSPC 


1992 


15-25/Apr 


22223 


600196A01 


49.3 


on 


6.14 ±0.38 


4 


PSPC 




17-27/Apr 


9443 


600195A01 


21.9 


on 


5.36 ±0.38 


5 


HRI 




23-24/Oct 


1096 


400237N00 


10.3 


off 


< 0.844 


6 


PSPC 




6/Dec 


3561 


6O0455NO0 


50.5 


off" 


- 


7 


PSPC 


1993 


29-30/Mar 


5214 


400300N00 


55.0 


ofF'''^ 


- 


8 


PSPC 




10-25/ Apr 


14207 


600452N00 


54.7 


Off'^ 


- 


9 


PSPC 




16-22/Apr 


1721 


600455A01 


50.5 


ofr"''= 


- 


10 


HRI 




17/Apr/93 


1167 


400237A01 


10.3 


off 


< 0.769 


11 


PSPC 




lO/May/93 


17593 


600453N00 


19.2 


off 


< 0.828 


12 


PSPC 




712-13/May 


4902 


500142N00 


52.6 


off 


< 3.91 


13 


PSPC 




1-14/Oct 


16663 


600452A01 


54.7 


on 


9.50 ± 0.70 


14 


PSPC 




1-9/Oct 


7199 


400300A01 


55.0 


on 


8.4 ±1.5 


15 


PSPC 




7-10/Oct 


4595 


600455A02 


50.5 


on 


16.4 ± 1.6 


16 


PSPC 




14-29/Oct 


20845 


500250N00 


53.5 


on 


9.95 ± 0.55 


17 


HRI 


1994 


1/Apr 


1301 


400237A02 


10.3 


off 


< 0.701 


18 


PSPC 




4-5/May 


4129 


600455A03 


50.5 


off 


< 0.395 


19 


PSPC 




5/May 


4111 


400300A02 


55.0 


off 


< 3.96 


20 


HRI 


1996 


26/Apr 


2006 


600810N00 


3.2 


off 


< 0.451 


21 


HRI 




10/June 


4720 


600810N00 


3.2 


on 


4.06 ± 0.32 



Table Notes: 

" The nearby source (RX J0054.9-7210) is active 

^ The nearby source (RX J0054. 9-7210) is active. The source state is derived from the nearby observation at smaller off-axis 
angle. 

" The source location is at the detector boundary, but no excess is seen 

^ The exposures listed are the nominal values. Close to the PSPC support structure the exposure is much reduced. The 
count rates are derived using an exposure map. 
^ Upper limits are given at 3 sigma 



Table 2. Optical magnitudes/colours of the two main optical candidates for IWGA J0053.8— 7226 on 1995 Nov 25/26 



Star 


U 


B 


V 


R 


I 


(U-B) 


(B-V) 


(V-R) 


(v-i) 


(R-I) 


Be star component of A 


14.004 


14.876 


14.927 


15.000 


14.995 


-0.87 


-0.05 


-0.07 


-0.07 


0.01 




(±0.011) 


(±0.016) 


(±0.023) 


(±0.024) 


(±0.024) 


(±0.03) 


(±0.04) 


(±0.05) 


(±0.05) 


(±0.05) 


F-star component of A 


14.553 


14.444 


14.075 


13.730 


13.412 


0.109 


0.369 


0.346 


0.664 


0.318 




(±0.016) 


(±0.007) 


(±0.006) 


(±0.005) 


(±0.005) 


(±0.023) 


(±0.013) 


(±0.011) 


(±0.011) 


(±0.010) 


Star B 


13.833 


14.637 


14.881 


14.725 


14.711 


-0.80 


-0.25 


-0.16 


0.17 


0.01 




(±0.016) 


(±0.007) 


(±0.006) 


(±0.005) 


(±0.005) 


(±0.023) 


(±0.013) 


(±0.011) 


(±0.011) 


(±0.010) 



Star A and indicate that they are separated by 2.1 arcsec- 
onds in an E-W direction. We give the photometric results 
in Table 2 for the stars labelled A and B in Figure 3. 

We find that Star B and one component of Star A arc 
very blue, and lie with a group of several other blue stars off 
the normal locus for stars in the two-colour {U—B)vs{B—V) 
diagram, implying some reddening. An E{B — V) of 0.15- 
0.20 would move the stars back onto the normal locus for 
early B-stars in the two colour diagram, namely a (B — V)o 
~-0.20 to -0.25. The fainter component of Star A has a 
{B — V) colour of 0.37 and is consistent with an early F 
spectral type, while the {V — R) and {V — I) colours imply 
a somewhat later type (F8). Our CCD spectra (Section 5) 
are more consistent with the earlier spectral type, so we 
adopt F2 for the this star's type. 



4 FOLLOWUP OPTICAL AND IR 
PHOTOMETRY 

Subsequent optical (BVRc) photometry was performed on 
1996 October 5/6. These CCD observations were also under- 
taken on the SAAO 1.0-m telescope, using the same (Tek8) 
camera, but employing the 3 x "Shara" focal reducer (kindly 
provided by Dr M. Shara, STScI), giving an effective field 
size of 9 X 9 arcmin, with a scale of 1.05 arcsec pixeP^ (this 
was utilised for a parallel program of optical identifications 
necessitating a larger field coverage). No [/-band observar 
tions were possible, due to the poor throughput of the focal 
reducer at this wavelength. Because of the plate scale, the 
close pair of objects labelled Star A in Figure 3 were not 
well resolved. In fact only for the B and V frames was it 
possible to derive PSF fits to both stellar images, although 
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Table 3. IR magnitudes for Star A and its nearby companion 
Star Date J H K 



CoTnbinsd^ 


Nov 95 


13.18 


12.93 


12.84 






(s.e. 0.05) 


(s.e. 0.05) 


(s.e. 0.07) 




Oct 96 


12.69 


12.47 


12.45 






(s.e. 0.10) 


(s.e. 0.07) 


(s.e. 0.05) 


Deconvolved: 










Nearby F-star 




13.36 


13.11 


13.08 


Be star 


Nov 95 


15.22 


14.97 


14.60 


Be star 


Oct 96 


13.53 


13.35 


13.34 



the results are far from reliable. The aperture magnitude of 
both stars (i.e. using an aperture large enough to include 
both) was V = 13.64, compared to a combined V — 13.667 
for 1995 November photometry, indicating very similar com- 
bined flux, and therefore little optical variability. The indi- 
vidually fltted PSF magnitudes were rather different, how- 
ever, no doubt a consequence of severe undersampling of the 
images, plus the fact that there are in fact 2 fainter stars 
within several arcseconds (Star C and F), i.e. four stars in 
total. Indeed our 1995 November observations resolved all 4 
stars. 

IR (JHK) photometry of Star A was performed on 
1995 November 14/15 and 1996 October 2/3 using the Mklll 
IR photometer on the SAAO 1.9-m telescope. A chopping 
secondary defines two effective apertures, one containing the 
source, the other background. After a short 10 sec integra- 
tion, the telescope "nods" such that the source is within the 
aperture previously defined as background. Thus two "off- 
source" positions are used, in order to minimise any system- 
atic errors. Observations are repeated, alternately nodding 
to the alternative background positions, until a sufficient 
S/N ratio is obtained. Most of the observations consisted 
of 5 observation modules of 40 s integrations, or a total of 
200 s. For both sets of observations, aperture sizes of 6 and 
9 arcsecond diameter were used, respectively, for the 1995 
and 1996 observations. 

The Ji?A' magnitudes are therefore of the combined light 
of both objects representing Star A (in fact also including 
the two other fainter, by 2-3 mags at V, neighbouring stars), 
and they are given in Table 3 with their standard internal 
errors. We were careful to observe a nearby JHK standard 
close in time (~40 min) to the IWGA J0053.8-7226 obser- 
vation. Therefore the true errors should not be too different 
to those quoted. Included are the deconvolved magnitudes 
for both components of Star A. These deconvolutions were 
made assuming that the IR colours, (V-K), (H-K) and (J- 
K), of the F-star could be determined on the basis of the ob- 
served (B-V), using the relationships between the SAAO IR 
JHK and Johnson-Cousin's UBVRcIc systems, determined 
by Caldwell et al. (1993). For the colours in question, these 
relationships are independent of luminosity class. From the 
observed V and predicted IR colours ( V-K, H-K, J-K) for 
the F-type contaminating star, we derived the JHK magni- 
tudes of both components of Star A. 

For the 1996 observations, the derived J-K for the Be 
star component of Star A is 0.19, and assuming E(J-K) 
~0.52 E(B-V), leads to a value oi {J - K)o ^ 0.11 (for E(B- 
V) — 0.15). This value can be compared to that expected 



for a late-O/early-B star (~-0.2), indicating that IWGA 
J0053. 8-7226 exhibits an IR excess of ~0.3 mags. Such an 
excess is a common feature of Be/X-ray binaries (e.g. Coe 
et al. 1994), and is probably related to the amount of mass 
present in the Be star's circumstellar disc. 



5 OPTICAL SPECTRA 

Our discovery spectrum, taken with the RPCS, lead to the 
identification of the Balmer emission (in Ha and weakly 
in H/3) in the eastern component of the close pair labelled 
Star A in Figure 3. However, the RPCS was unable to ob- 
tain spectra of both components separately. This was sub- 
sequently achieved using the SITe2 CCD detector on the 
Cassegrain spectrograph, with the same grating and wave- 
length coverage. The observations (2 x 600 s) were kindly 
undertaken by H. Winkler on 1997 October 2/3, in photo- 
metric conditions. The flux standard star EG 21 was ob- 
served immediately after the IWGA J0053. 8-7226 observa- 
tions, and the data reduced using IRAF[|| tasks in ccdproc, 
onedspec and twodspec. 

The close proximity of the two objects in Star A meant 
that the individual spectra were hardly resolved in the spa- 
tial direction (i.e. perpendicular to the slit) on the CCD, 
having a typical cross section of FWHM ~3 pixels. Using the 
apall program, we experimented extracting the spectrum 
using a range of apertures. The Hq emission was clearly 
visible over several CCD rows, which we used as an initial 
indicator for which rows to extract for the Be star spectrum. 
The best results were achieved by optimally extracting the 
spectra over a range of rows, starting 6 pixels to one side of 
the peak in the cross section, to within 1 pixel of this peak 
(or 5 pixels in total) . Likewise, the contaminating nearby F- 
star was extracted over the remaining pixels in the proflle, 
i.e. 6 pixels on the other side of the peak to the limit of the 
Be star extraction (7 pixels). The spectra were flrst summed 
over 30 pixels at a time in the dispersion direction, to obtain 
the centroid (in the spatial direction) of the spectrum, and 
then traced using a cubic spline. 

Observations of Star B were subsequently carried out 
using the same telescope and instrument set up on 19/20 
August 1998. 

Our results are shown in Fig. 4 where the spectra of 
both Stars A and B are shown. The spectrum of the F star 
near to Star A is shown in Figure 5. Although sky back- 
ground has been subtracted, there has been no attempt to 
remove the telluric lines (e.g. the A-band). The nearby F- 
star shows strong Ca H & K and Balmer absorption lines, 
on a blue continuum. Judging from the relative strength of 
the G-band with respect to H7, an early F spectral type is 
implied, consistent with the observed {B — V) of 0.37. The 
Ha absorption line maybe somewhat flUed in from contami- 
nation by the Be component to Star A. Indeed both spectra 
are likely to suffer from mutual contamination at some level, 
although the results are far superior to the RPCS spectra. 
We fltted Gaussian profiles to the emission lines for Stars 
A and B in both the RPCS discovery spectrum and the 



t IRAF is distributed by NOAO, which is operated by AURA, 
Inc., under contract to the NSF. 
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Figure 4. Top panel: spectrum of the eastern component of the 
star labelled as A in Fig. 3. Lower panel: spectrum of Star B. 




" 4000 4600 5000 6500 6000 6500 7000 7500 
Wavelength (A) 



Figure 5. Spectrum of the F2 star which is the western compo- 
nent of the double labelled as A in Fig 3. 

CCD spectra, and the results are summarised in Table 4. 
The radial velocity values arc consistent with membership 
of the SMC, which has a mean systemic value of 166 ± 3 
km s~^ (Feast 1961). The same is true for the nearby F2 star, 
which has a similar velocity shifted Balmer absorption lines. 

The remaining stars observed in, or near, the error circle 
are very unlikely to be the counterpart to the X-ray source. 
Stars C and E are clear A type stars. Star D is an F type. 
Star F a carbon star and G is probably a late K type, proba- 
bly dwarf star (and also a significant distance from the X-ray 
position). 



6 DISCUSSION 

Apart from Ha, H/3 and the other Balmer lines (seen in 
absorption H7, HJ, He, H8), there appears to be no other 
significant features in the spectra of Stars A and B. Regret- 
tably the S /N ratio of the spectrum is insufficient to obtain 
further spectral typing information. 

However, for the Be component of Star A we can use 
the contemporaneous optical and IR photometry from 1995 
November to model the spectral shape. Assuming a value 




WavBlength (A) / 10* 



Figure 6. Photometric data points together with the best-fit Ku- 
rucz model for a B1-B2 III-V star (Teff = 22000 ± 2000K; logs = 
4.0). 

of E(B V) = 0.15, we de-reddened our data and fitted a 
Kurucz model atmosphere to the photometry. Normalizing 
the model to our U band point results in an acceptable fit 
for T = 22000 ±2000 K and \ogg = 4.0. This fit is illustrated 
in Fig. 6, from which the IR excess over the stellar model 
is clearly visible. This temperature range corresponds to a 
B1-B2 spectral type. 

Unfortunately we have no IR data for Star B to demon- 
strate that the IR excess for that object also exists. How- 
ever, the (U-B) and (B V) optical colours are very similar 
to those of Star B indicating a similar spectral type. 

In the case of the Be star component to Star A mem- 
bership of the SMC is clearly established. Hence we can 
estimate the absolute magnitude of this system using the 
observed apparent V magnitude in Table 2, correcting for 
an assumed absorption of E(B-V) ~0.2, implying an Ay 
^0.6. This gives my ~14.3, and for a distance modulus of 
{mv - M„) = 18.9 (Westerlund 1997), implies My = -4.6. 
This suggests a luminosity class III-V. For the nearby con- 
taminating F2 star, the absolute magnitude is estimate to 
be Mv,c = -4.8 to -5.4, respectively, for an E{B — V) 
in the range 0.0 to 0.2. This star is therefore most likely a 
supergiant F-star. 

The degree of reddening seems consistent with the col- 
umn density, Nh, determined from the ROSAT observa- 
tions. An E{B — V) of 0.15 implies a column of Nh = 
1.4 X lO^^atoms cm~'^, which compares favourably with 
the value obtained from the X-ray spectrum (Section 2), 
namely 2.3± 1.1 x 10^^ atoms cm~^. These column densities 
are all much larger than the column density to the SMC 
(6 X lO^^cm"^), implying a local source of reddening consis- 
tent with a circumstellar disk around the Be star. 

In Section 2 we estimated the unabsorbed 0.1-2.0 keV 
X-ray flux as ~7.5 xlO"^"^ ergs cm^^ s^^. Using the same 
power- law model, we calculate the unabsorbed 2-10 keV flux 
as ~3.2 xl0~^^ ergs cm~^ s~^, which we use to determine 
the luminosity ratio, Lx{2 - WkeV) : Lopt (3000 - 7000A), 
the ratio used by Bradt & McClintock (1983) in their sur- 
vey of X-ray source optical counterparts. Our estimate of 
fopt,unabs ~2 xlO"^'^ crgs cm~^ s~^ for Star A comes from 
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Table 4. Gaussian fits to Ha and H/3 



Object Date Radial Velocity Flux E.W. FWHM 

(km s-l ) ( ergs cm-^ s"!) (A) (A) 



Star A Ha: 

(Be star 22/23 Sep 95 177 3.6 x10"" -13.1 8.7 

component) 2/3 Oct 97 157 3.7 x10"" -11.6 9.2 

(±15) (±0.3 xlO"") (±0.9) (±0.7) 

H/3: 

2/3 Oct 97 186 4.1 xl0"i5 -0.8 5.4 

(±50) (±1.2 Xl0"i6 (±0.2) (±2.1) 



Star B 



Ha: 

19/20 Aug 98 
H/3: 

19/20 Aug 98 



5.8 xlQ-" 
(±0.3 xlO"") 

1.3 xlO"" 
(±1.2 xl0"i5 



-15.3 
(±0.9) 

-1.76 
(±0.2) 



10.7 
(±0.7) 

6.9 
(±2.1) 



the optical spectrum (Fig. 4), where wc integrated the flux 
over the interval 4000- 7000 A, and extrapolated linearly to 
3000A. The derived Lx/Lopt (= fx,unabs/ fopt,unabs) ratio of 
^^0.01-0.02 is very typical of HMXRBs, and the more active 
Be/X-ray binaries. The ratio is some 4 orders of magnitude 
greater than that expected from pure coronal emission from 
early type stars (Caillault & Helfand 1985), and consistent 
with accretion powered X-ray emission. The inferred lumi- 
nosity of IWGA J0053.8-7226 , Lx (2-lOkcV) = 1.6x10^^ 
ergs s^^, is lower than X-ray "high-state" luminosities often 
seen in transient Be/X-ray systems, which are typically in 
the range ergs s . The value is far more typical 

of the quiescent luminosity produced by stellar wind accre- 
tion (Davidson & Ostrikcr 1973). 

It is interesting to note that IWGA J0053.8-7226 does 
not appear in the ROSAT HRI survey of the SMC (Gowley 
et al. 1997). An upper limit of 0.001 cts s"'^ was placed on 
the flux, whereas it was expected to bo 0.007 cts s"^ based 
on the Einstein count rate reported in Wang and Wu (1992). 
This clearly establishes the source as variable, since it was 
observed by both Einstein and the ROSAT PSPG (during 
the RASS). 

IWGA J0053.8 7226 is clearly a transient source, with 
excellent bona fides for a Be/neutron star system. The X- 
ray, optical and IR variability (at least for one candidate), 
combined with the 92 s X-ray pulse period (Corbet et al. 
1997), convincingly demonstrates that IWGA J0053.8 7226 
is a Be/neutron star binary system, though we are unable 
to distinguish between the two possible optical counterparts 
at this stage. Only if correlated optical/X-ray variability is 
seen, or a much more accurate X-ray position determined, 
will this point bo resolved. This makes it the S*** such sys- 
tem discovered in the SMC, following the discovery of RX 
J0117.6-7330 (Clark, RemUlard and Woo 1997 and Coe et 
al 1998). 

We also note that a signiflcant number of X-ray pulsars 
have recently been discovered in this region of the SMC. 
In addition to the source RX J0051.8-7231, from which 8.9- 
second pulsations were reported by Israel et al. (1997), the 



sources AX J0051-722 (91 second pulsations, Corbet et al. 
1998), IWGA J0054.9-7226 (59 second pulsations, Marshall 
et al. 1998) and RX J0049. 1-7250 (75 second pulsations, 
Yokogawa and Koyama 1998) have all been discovered in 
recently within a 40 arcminute region of the bar of the SMC. 
We have identified the optical counterparts to these X-ray 
pulsars to be SMC Be stars through similar optical programs 
as described here (Stevens, Coe and Buckley, 1999). 

It is likely that the number of recent discoveries is sim- 
ply due to the intensity of observations of this region follow- 
ing the discovery of the first pulsar. Future observations of 
the whole of the SMC with the new generation of X-ray tele- 
scoi)os should reveal whether this region really does contain 
an anomalous high density of X-ray pulsars. 
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